An Onsager-like Relation for the Lattice Boltzmann Method 
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An Onsager-like relation is proposed as a new criterion for constructing and analysing the lattice 
Boltzmann (LB) method. For LB models obeying the relation, we analyse their linearized stability, 
establish their diffusive limit, and find new constraints for those with free parameters. The new 
relation seems of fundamental importance for the LB method. 
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The lattice Boltzmann (LB) method is an effective and 
viable tool for simulating complex fluid flow problems. 
Historically, the method was evolved from the lattice gas 
automata [1, . It is now well understood as a special 
discretization of the Boltzmann equation The ki- 
netic origin enables the method to naturally accomodate 
a variety of boundary conditions for flows with complex 
geometry. Besides its effectiveness and simplicity (see 
Eq. ^ below), the method has a clear and solid physi- 
cal interpretation [3] ■ These advantages make it have 
vast applications. Indeed, the method has been success- 
fully used in a wide spectrum of areas including turbu- 
lent flows, microflows, multi-phase and multi-component 
flows, and particulate suspensions. It is becoming a seri- 
ous alternative to traditional computational methods in 
fluid dynamics. We refer to [3, [3, for reviews of the 
method and its applications. 

In spite of the vast and successful applications, the nu- 
merical stability of the method has not been well under- 
stood but has attracted much attention [3, 0, @, [13, 12, 



15l . Il8l . l2d . |21| . In this work, we intend to understand the 
LB method from the viewpoint of nonequilibrium ther- 
modynamics [llj . This understanding will provide a 
new insight into the method and a novel approach to 
study the stability. 

The general form of the LB method is 

/,(x^+c„< + l) = /,(x^,<) + J,(f(x^,<)) (1) 

for i = 1,2,--- and fi in a finite or countable set. 
Here fi — fi{'x.^,t) is the probability of finding a fluid 
particle, at site x^ G and time i > 0, that travels 
with velocity S W^, Ji{^) is the z-th collision term, and 
^ — ifii f2, ■ ■ ■ , Jn)'^ with the superscript T denoting the 
transpose operation. We will often write 

J(f) = (Jl(f),J2(f),-- - ,JNii)f. 

Eq. |T]) is fixed by specifying (or constructing) the 
discrete velocities and collision terms Ji(f). This is 
usually done by requiring ([1]) to possess the following 
properties: consistency with certain conservation laws, 
Galilean invariance, isotropy, small or no compressible 
effects, velocity-independent pressure, and so on. See 
for details. 



In this letter, we will propose a new requirement for 
the construction. This requirement has its origin in non- 
equilibrium thermodynamics. According to de Groot and 
Mazur [ol, nonequilibrium thermodynamics classifies ir- 
reversible phenomena into three types: scalar, vectorial 
and tensorial processes. Typical examples are chemical 
reaction, heat conduction and viscous fiows, respectively. 

Our starting point is to understand that the A'' equa- 
tions in ll]) describe N scalar irreversible processes and 
the collision terms Ji(f) are the corresponding irre- 
versible fiuxes. Moreover, the thermodynamic forces are 
identified as the gradient of a possibly-existing entropy- 
density function H{f) [l^. With such identifications, 
the Onsager reciprocal relation reads as 



J(f) = 5(f)i//(f), 5(f) = 5(f)^. 
The null space of matrix ^(f ) is independent of f . 



(2) 



Here Hf{f) stands for the gradient oi H = H{i) with 
respect to f , that is, 

Hf{t) := (%,%,... ,(}/„)^iJ(f). 

Let us mention that the gradient of the entropy-density 
function H = H{f) has been used in to construct 

entropic LB methods. 

Remark that, in the classical Onsager relation 0, the 
syinmetric matrix S = S{f) is constant. Recently in 
the author looses the constancy requirement and 
proposes (2) for general scalar processes. The new inde- 
pendence in (2) expresses the fact that physical laws of 
conservation hold true, no matter what state the under- 
lying thermodynamical system is in (equilibrium, non- 
equilibrium, and so on). In addition, the second law of 
thermodynamics requires that the symmetric matrix ^(f ) 
have a sign (non-positive). 

If an LB model admits the relation ([2|), one can eas- 
ily show that it allows an H theorem with H{i) as the 
entropy-density function (see [13, [3]). Therefore, the 
relation ^ is not realistic, for many widely used LB 
models were shown in 21 1 not to admit an H the- 
orem. On the other hand, it follows from (2) that 
S{QHf{f.,) = J{f.,) = for f* satisfying J{f.,) = (in 
equilibrium). Thus, Hf{i^) is in the null space of S'(f*) 
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and thereby in that of 5'(f), since the null space of 5(f) 
is independent of f . Therefore, we have 

S{i)Hf{U) = 

for any f . Now we differentiate the two sides of Eq. ^ 
with respect to f and compute at f* to obtain 



d_ 

Of 



[Si{)Hf{{. 



|f=f. 



(3) 



Here Hff{i) is the Hessian of the entropy function H{i). 

From Lemma 2.1 in j22| we know that, if the discrete- 
velocity set is such that ^ Cj whenever i ^ j, the 
possibly-existing entropy-density functions H{i) for ((!]) 
must be of the following form 

N 



1=1 

That is, H{f) does not contain any cross-terms. Thus, 
the Hessian Hff{t) must be a diagonal and positive- 
definite matrix. The positive-definiteness follows from 
the strict convexity of H(f) as an entropy-densitiy func- 
tion. 

Inspired by these considerations, we propose the fol- 
lowing notion. The LB method ([1]) is said to obey an 
Onsager-like relation at f = f, satisfying J(f,) = 0, if 
the Jacobian matrix of J(f ) evaluated at f = f* can be 
decomposed as 



(4) 



where ^(f,) is a symmetric non-positive-definite matrix 
and DJf^ is a diagonal positive-definite matrix. 

In [3, |8| , we verified that many LB models admit the 
following structure: there is an invertible N x A^-matrix 
P such that P^P is diagonal and 

FJ/(f,) = -diag(Ai,A2,--- ,Ajv)P (5) 

with the Ai's non-negative. In this situation, we have 

J/(f,) = -F-Miag(Ai, A2, • • • , Xn)P- 



This is the relation ([J]), because P^ P is diagonal and 
P~^diag(Ai, A2, • • • , Xn)P~^ is symmetric. 

It is interesting to note that the relation ^ is equiva- 
lent to dni). To see this, we set A = ^D(f,). As AS'(f,)A 
is symmetric and non-positive definite, there is an or- 
thogonal matrix U such that 

A5(f,)A = -C/^diag(Ai,A2,-- - ,Xn)U 

and the A^'s are non-negative. Thus, it follows from Eq. 
dH) that 

Jf{U) = S{U)A^ = -A-i[/^diag(Ai, A2, • • ■ , \n)UK. 



Set P = UK and notice that U'^ = [/^. Then 
pTp = a2 = D{U) is diagonal and PJ/(f,) = 
— diag(Ai, A2, • ■ ■ ,Xn)P- Hence the relation ([4]) also im- 
plies ©. 

Generally speaking, the invertible matrix P in ([5]) dif- 
fers from the transformation matrix M used in construct- 
ing the multiple-relaxation-time LB method [i, 12|. The 
latter is such that MM'^ is diagonal. 

In the rest of this article, we present some consequences 
of the Onsager-like relation ([4]) or its equivalent version 
®. 

First, the relation (|4]) can be used as an analytic cri- 
terion to fix LB models with free parameters. In [3], we 
require that a number of existing LB models with free 
parameters obey the equivalent relation ([5|) and find new 
constraints for the free parameters. In this way, either 
the free parameters are fixed or their freedom is reduced 
considerably. The results coincide with those obtained 
through numerical tests and/or guesswork. 

An example is the D2Q9 model with two parameters 
a and (3 [15|]. Here d = 2, = 9, the discrete velocities 
are Cg = 0, 

{c. :l<i<4}={(±l,0f,(0,±l)^}, 
{c,:5<i<8}-{(±l,±l)^}, 



and the coUsion terms are 



with r a relaxation time. Furthermore, 
o 

Jg (ri, v) =an - -n\^^\ , 



for i = 1,2,3,4, 

/f«(n, v) = ^ —n H vc^ -I- -(vCiV vr 

Jz \ ^ ) 4 12 8 24' 

for i = 5, 6, 7, 8, and 

9 9 
n = ^fz, nv = ^c,/i. 

i=l i=l 

In we verified that the above parametrized D2Q9 
model obeys the equivalent version ^ if 



a e (0, 1) and a + 5/3=1. 



(6) 



This is satisfied by both (a,/3) = (4/9,1/9) used in the 
well-known D2Q9 model jj] and (a, l3) = (2/7, 1/7) used 
in Moreover, the numerical tests in |2| support the 
other choices of (a,/3) satisfying condition ([6]). 

Next, we show that the equivalent version (O provides 
a convenient setting to analyse the linearized stability of 
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the LB method ([T]) . Let f* be a uniform equihbrium state 
and f the fluctuation. The hnearized LB method is 



f^i+l) = f(x^,i) + J/(f,)f(x^,t), 



(7) 



where 



fp(i) = (/l(x^+Cl,t),-- - ,/Ar(Xp +Cjv,t))^. 

For the sake of simphcity, we consider only the periodic 
initial data for Q, where fi ranges in a finite set. For 
other cases, see [9|. 

Assume the Onsager-like relation dU holds at the equi- 
librium state f.^. Then we have the equivalent relation 
([5]) . Multiplying the linearized equation ([7]) with P from 
the left gives 

Pif,{t + 1) = diag(l - Ai, 1 - A2, • • • , 1 - A7v)Pf (x^, i)- 
Thus, if 

A. G [0, 2] V i, (8) 

it is obvious that 

\PUt + l)\' <\Pi{x^,t)\\ (9) 

Here |f| is the Euclidean length of the iV-vector f. On 
the other hand, since P^P is a diagonal matrix, say 
diag(ai, 02, • • ■ , a^r), it follows that 



\Pi^{t + =i^{t + ifP'^Pi^it + 1) 

N 



Because x„ + is a lattice node, we deduce that 



TV 



N 



1=1 



= ^|Pf(x^,i + l)|^ 



method. Moreover, the derivation is different from those 
in [lo| . does not involve any entropy-density function, 
and works for LB models violating H theorems. Finally, 
our approach can be easily extended to other cases, even 
with boundaries 

Finally, we mention that the Onsager-like relation ^ 
ensures the validity of the formal diffusive limit for the 
continuous version of the LB model ([1]) parametrized 
with e > 0: 



^tf^ 



1 



1 



Mf). 



(11) 



It is well-known (see, e.g., [8|, |15[) that this system 
of parametrized partial differential equations is related 
closely to the LB method (H]). Suppose the discrete- 
velocity set is symmetric in the sense that 

{c, : I = 1, 2, • • • , iV} = {-c, : ^ = 1, 2, • • • , AT}, 

the parametrized model (|lip is consistent with the incom- 
pressible Navier-Stokes equaton and obeys the equivalent 
version ([5]) at a quiescent state f*. We rigorously veri- 
fied in 8] that the solution /f (x, t) to (fTTj) with prepared 
initial data satisfies 



N 



i=l 



//(x,t)c, = ev(x,t)-t-0(e3 



(12) 



as e goes to zero. Here v(x, t) is the velocity of the fluid 
under consideration. The details are given in Q. 

In conclusion, we have proposed an Onsager-like rela- 
tion ^ and its equivalent version ([5]) as a new require- 
ment for analysing and constructing LB models. For 
LB models obeying the Onsager-like relation, we anal- 
yse their linearized stability, established their diffusive 
limit, and found new constraints for those with free pa- 
rameters. On the basis of these consequences, we believe 
that the new relation Q is of fundamental importance 
for the LB method. Finally, we expect that the Onsager- 
like relation ^ can be used as a guide to construct stable 
LB models. 



Now, summing up Eq. ^ over all ^ (in the flnite set!) 
gives the following inequality 



^|Pf(a;^,t+l)p<^|Pf(x^,t)| 



(10) 



This simply means the stability of the linearized LB 
method ([7]) under the condition ([8]). Note that for the 
D2Q9 model above, ([8]) is nothing but the well-known 
condition r > 1/2. 

Our derivation of the inequality (jlOp does not involve 
the von Neumann stability analysis, which was used in 
the previous works [l^, [3, [l3| on the stability of the LB 
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